Influence of Coulomb and Phonon Interaction on the Exciton Formation Dynamics in 

Semiconductor Heterostructures 



Walter Hoyer,^ Mackillo Kira/ and Stcphan W. Koch^ 

'Department of Physics and Material Sciences Center, 
Philipps University, Renthof 5, D-35032 Marburg, Germani]^ 
(Dated: February 1, 2008) 

' A microscopic theory is developed to analyze the dynamics of exciton formation out of incoherent 

' carriers in semiconductor heterostructures. The carrier Coulomb and phonon interaction is included 

^ ^ , consistently. A cluster expansion method is used to systematically truncate the hierarchy problem. 

^ ■ By including all correlations up to the four-point (i.e. two-particle) level, the fundamental fermionic 

O ' substructure of excitons is fully included. The analysis shows that the exciton formation is an 

I intricate process where Coulomb correlations rapidly build up on a picosecond time scale while 

phonon dynamics leads to true exciton formation on a slow nanosecond time scale. 

CN ' PACS numbers: 71.35.-y,71.10.-w,73.21.Hb 

Keywords: exciton formation, many-body correlations 

^ ■ 

ill ! I. INTRODUCTION 



I 



The linear absorption spectrum of an intrinsic direct-gap semiconductor in its ground state shows excitonic reso- 
nances energetically below the fundamental band-gap energy. These resonances in the linear optical polarization are 
a consequence of the attractive Coulomb interaction between a conduction-band electron and a valence-band hole. 
Mathematically, the eigenvalue problem of an exciton iruthe so-called Wannier limit is identical to that of the hydrogen 
atom, such that all relevant properties are well known. EJ 
Q ■ However, the appearance of excitonic signatures in optical spectra does not necessarily imply the population of 
O ', excitonic states. Even if the optical excitation is directly resonant with the exciton energy, the optical excitation 
always induces optical polarization in the semiconductor. This polarization, respectively the absolute square of the 
I ' polarization, is sometimes related to coherent populations; however, it does not yield information about the truly 
^ , incoherent population of single particle or pair states. 
^sD ' Furthermore, the optically induced polarization and the other coherences decay away typically orua picosecond 
OO ' time scale|-ji|ie to: i) excitation-induced dephasing rosttlting from the carrier-carrier Coulomb scatteringpa, ii) phonon 



' scatteringEroO, and iii) the finite radiative lifetimda'QQ in confined semiconductor structures. However, after the 
' I [ disappearance of the optical polarization the material is typically not back in its ground state since carrier densities 
■ and incoherent correlations remain in the system and continue their many-body dynamics for several nanoseconds. The 
life time of these incoherent excitations is limited only by slow radiative and non-radiative recombination processes. In 
-j—; ideal semiconductors, only spontaneous recombination determines the ultimate lifetime of the excited charge carriers. 

In this context, an old but still open question concerns the nature of the incoherent Coulomb correlated populations, 
^ , in particular the conditions under which a significant part of the clcctron-holc excitations exist in the form of bound 
I ■ pairs, i.e. Wannier excitons. Although it is not even clear a priori how to count bound states since a rigorous exciton 
■ O number operator does not existcl, one wants to understand the quantum statistical prcmacties of incoherent excitons, 
^ i their distribution function, and possible bosonic as well as Bose condensation aspects.t3Ej 

^ ' So far, several papers disciias-Uie relaxation dynamics of excitons due to phonons together with a direct bosonic 
• • : approximation for the excitorO'tj. In other approaches, exciton formation after coherent excitation has been investi- 
. ^ gated, howevein-|pstricting the treatment to weak excitation conditions due to a low-order expansion in terms of the 
K> , exciting pulse.tllE-3 Since these approaches do not involve fermionic carrier densities and carrier-carrier correlations, 
5^ ' they cannot resolve how much of the excitation is transferred directly into those correlations and densities. Typically, 
] in bosonic models, exciton contributions are separated into bound and continuum excitons. In this description, one 
cannot determine the plasma carrier density since the continuum excitons determined from the two-particle reduced 
density matrix are not trivially related to the carrier densities given by the one-particle reduced density matrix. 

In this paper, we present and utilize a microscopic theory for the description of the totally incoherent regime. In this 
approach, fermionic carriers, incoherent exciton populations, carrier-carrier correlations, as well as phonon-assisted 
correlations are treated at the same level. As a result, we are able to evaluate the excitonic correlations, formation 
rates, distribution functions, etc. for different temperatures and carrier densities without neglecting the underlying 
fermionic anti-symmetry. 

After presenting the t heor etical background and the full equations in the incoherent limit in Sec. ||, we discuss our 



numerical results in Sec. III. This is done by viewing the build-up of correlations and monitoring the growth of the 
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electron-hole pair-correlation function. We investigate the issue of energy conservation and derive an approximate 
adiabatic solution which offers the possibility to obtain an intuitive interpretation of the exciton formation process. 
The results are summarized in Sec. |^ and some technical details of the calculations are presented in the appendices. 



II. INCOHERENT EXCITON AND PLASMA DYNAMICS 

In contrast to earlier investigations of exciton formation after resonant excitationElO, we focus in this paper on 
an initially completely incoherent configuration. Such a situation can be realized experimentally either by pulsed 
continuum excitation, after the initial coherences have decayed, of with a current injection of carriers since this 
process does not induce optical coherences. In all numerical investigations in this paper, we therefore start with 
vanishing correlations and finite carrier densities and view the evolution of the system developing from that initial 
condition. 



A. The Hamiltonian 



The starting point of our theoretical description is the full Hamiltonian which describes the free motion of non- 
interacting carriers and phonons as well as the interaction between those quasi-particles. In this sectioB.|-we briefly 
outline the Hamiltonian which is used throughout the paper; for more detailed derivations, see e.g. Refs.lllll3. 

In semiconductor heterostructures, the microscopic properties of the carriers can be described with the fermionic field 
operator a\,k {a\ f.) annihilating (creating) an electron in band A with wave vector k along the heterostructure. Here, 
in general, the band index A may include different bands, subbands, and spins. For sufficiently narrow confinement, 
the carrier dynamics is restripted practically completely to the lowest subband and the Bloch functions follow from 
the envelope approximation.EJ The corresponding free Hamiltonian is 

iJkin = (fffe al kQck + el alj^ay^k^ , (1) 

k 

which includes one valence and one conduction band. The generalization to a multi-band system or the inclusion of spin 
is straight forward.BEll In general, the eigenenergies e'f/^ must be obtained from precise band structure calculations. 
For the investigation of near-bandgap optical features, we may use = e^. -I- Eq = |^ + Eq and = — ej! = — f;;^ 
where and m'* are the effective electron and hole masses. The unrenormalized bandgap is denoted with Eq^^ 
In order to include the Coulomb interaction between the electrons in different bands, we use the Hamiltoniaru'EJ 

iJc = ^ Y ^1 {^l,k(^l,k''^c,k'+qac,k~q + al ,^al i^,av.k'+qav.k-q + 2 al ,^alj^,ay,k'+qac.k-q^ , (2) 

k,k' ,q^0 

where Vq is the quantum-well or quantum-wire Coulomb matrix element. The first two terms in Eq. (H) lead to 
repulsive interaction between electrons within the same band whereas the last term gives rise to the attractive 
interaction between electrons and holes (missing valence band electrons) in different bands. 

In semiconductors, the carriers are additionally coupled to lattice vibrations (phonons). The noninteracting part 
of the phonon Hamiltonian is 

i/phon = 5]?if^|p| fel)p + ^) , (3) 



p 

where the bosonic operator Dp {D'^) annihilates (creates) a phonon in the state p. Here, p = {p,p±,a) labels the 
wave vector along and perpendicular to the heterostructure as well as the phonon branch. Corresponding to the three 
independent modes of sound waves in a solid, three branches of acoustic phonons always exist in a three-dimensional 
semiconductor. tj In lattices with more than one atom within a unit cell, additional optical phonon branches are 
present. 

The carrier-phonon interaction Hamiltonian can be expressed in the form 

Hp= Y ^lp±{^P,P±^l,k"-'^^k+p + Gl p^al ,^ay,k+p) + h.c. , (4) 

p,p±,k 

where Gp is the effective coupling matrix-element between a phonon and an electron in band A. Only the momentum 
component p along the heterostructure is conserved because the Bloch electrons are confined in perpendicular direction 
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p± . The form of Eq. (Q) reflects the fact that phonon interaction takes place within one band because phonons cannot 
provide the energy for an interband transition. The exact form of the interaction matrix element Gp[p\ depends on 
the interaction type which is included (i.e. deformation potential, piezo-electric coupling, etc.). For the treatment of 
acoustic phonons in the long- wavelength limit, the matrix element is given by a product between the plain matrix 

element Cp = ^ 2V^clap ^' ^^i^re p denotes the mass density and V the volumne of the semiconductor material, cla 
is the velocity of sound and D the deformation constant, and the form factor fR{p±) which is a consequence of the 
confinement of the electrons in the direction perpendicular to the semiconductor structure. By choosing parabolic 

confinement, we obtain = e 2 for a one-dimensional quantum wire characterized by the confinement 

potential 

The starting point of all our further investigations is the total Hamiltonian of Eq. 

Htot = ffkin + He + Hphon + Hp, (5) 

where both interaction parts of the total Hamiltonian lead to non-trivial coupling to higher order correlations. 



B. Hierarchy Problem 

Starting from Eq. (|5|), one can use the Heisenberg equation to compute the equations of motion for all relevant 
expectation values of interest. The simplest examples of incoherent one-particle expectation values are the two-point 
quantities describing the electron and hole distributions 

/l = (<fe«c,fe), = K-feflLi:) • (6) 

As is well known, the dynamics of does not yield a closed set of equations since both the Coulomb and phonon 
interaction Hamiltonians lead to coupling of the two-poititjauantities to either four-point quantities or mixed carrier- 
phonon operators. Schematically, the hierarchy problcmtltEl is described by 

^h^^{N)=Tm]+V[{N + l)] (7) 

where an iV-particle (i.e. 2Af-point) operator is coupled to (Af-l-l)-particle operators due to the many-body interaction. 
To solve this hierarchy problem rigorously, we apply a method known from quantum chemistry where the tru ncati i 
problem for many-electron wave functions has successfully been approached with the so-called cluster expansion 
There, the electronic wave function is divided into classes where electrons are: i) independent single particles (singlets), 
ii) coupled in pairs (doublets), iii) coupled in triplets, and iv) coupled in higher order clusters. The TV-particle wave 
function is then approximated by a suitable amount of coupled clusters by including the correct antisymmetry of the 
fermions. 

In semiconductors, the system properties can be evaluated from 2iV-point expectation values (A^) = 
{a\ ^ ..a\ ^ a^„^p„ ..Ojyj^^pj ), which determine the reduced density matrix in the Bloch basis. If the system con- 
tains exactly M particles, (M) fully describes the system properties. In this case, 2iV-point expectation values exist 
for all N < j\f. The simplest approximation scheme is provided by the Hartree-Fock approximation which implicitly 
assumes that the many-body system is described by a so-called Slater determinant of J\f independent single-particle 
wave functions. In this case, 2A'^-point expectation values can be expressed in terms of two-point expectation values. 
For example, a four-point expectation value in singlet (i.e. Hartree-Fock) approximation is given by 

^'^Ai,fei'^l2,fc2'^''2,P2'^I^l,Pl )S = (<^Ai,fci'^''l,Pl)('^A2,fc2'^''2,P2) ~ ('^Ai ,fei '^!^2 ,P2 ) ("^Aa ,fc2 '^'^l .Pi ) ' 

In the spirit of the original cluster expansion the system is thus fully described by uncorrelated singlets where each 
carrier behaves effectively like a single particle influenced by the meanfleld of all other particles. Higher order 
contributions are deflned recursively. For example, true two-particle correlations are obtained from 

^^'^Ai,fei'^A2,fc2'^'^2.P2'^I^l,Pl ) = ('^Ai,fci'^A2,fe2'^'^2,P2'3^!^l,Pl ) ~ ('^ Ai ,fci A2 ,^2 '^'^2 ,P2 '^^I'l ,Pl ) S , (9) 

i.e., by subtracting the singlet contribution from the full expectation value. The advantage of this factorization is 
its direct physical interpretation. By subtracjtiHg the single-particle contribution, the resulting correlated part really 
describes the true two-particle correlations fiSEj In the same way, truly correlated triplets and higher clusters can be 
defined by subtracting all lower-level contributions.E£l 
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If we only include plasma and pair-correlation effects in the analysis, an arbitrary iV-particle expectation value 
can be expressed consistently with the singlet-doublet approximation (N) « (iV)gD. This determines uniquely how 
the truncation of the hierarchy has to be performed; we only need to — and are allowed to — solve dynamics of 
(1) and A(2) because any arbitrary (A^)sd consists only of two-point expectation values and four-point correlations. 
According to Eq. (|^, we obtain 

= T,[{l)] + V^ms]+V^[A{2)],, (10) 
i?i^A(2) = r2[A(2}]+F2[(3)sD], (11) 

where 2^1(2) and V\(2) are known functionals defined by the specific form of the Heisenberg equation of motion. The 
consistent singlet-doublet-approximation is obtained when (3) is approximated by (3)sd; as a result, the infinite hier- 
archy is systematically truncated and Eqs. (p^)-(|ll|) are closed. In order to study exciton formation, we consequently 
have to know the exciton and carrier-carrier correlations 



= A(a^a^,ac,fc'+qa„,fe-,), (12) 
^g.k'.k ^ A(a^a^,a,,fe.+,a,,fe_,), (13) 

•h'' = '^{0'ik°'ik''^v,k'+qav,k-q), (14) 



in addition to the single-particle densities, Eq. (g). 

In geBetal, the {N)sd truncation fully includes the x^^^ limit of the so-called dynamics controlled truncation 
scheme. cllc3 Obviously, a similar relation holds as the number of clusters is increased; i.e. is a subset of the 

A^-particle cluster expansion. Furthermore, the cluster expansion can be used in regimes where x^''^^ methods fail, 
e.g., when the exciting field is strong or the system becomes fully inppherent. When (3) is approximated by including 
only the singlet part (3)s in Eq. (pi]) , also the second Born approachd is found to be a subset of the cluster expansion. 
Since in this case the doublet part is included only partially, the second Born scheme does not allow formation of 
bound states but rather describes microscopic scattering. This ideology can be directly generalized in order to include 
the triplets in the scattering level; for such an approach, the A (3) dynamics is solved but the resulting (4)sdt is 
approximated by (4)sd. This simplifies the numerical complexity of the triplet terms considerably and allows their 
analytic evaluation on the scattering level. 

The cluster expansion can also be used directly to classify and truncate carrier-phonon correlations using the formal 
equivalence between a boson operator and a pair of fermion operators. The operator equation of motion. 



= nQp^p^ Dp^p^ + g^p^a\ f^ax,k+q, (15) 

fl'phon+-f^P A,fc 



shows that one phonon operator is formally equivalent to a product of two electronic operators and that with every 
phonon absorption or emission process, an electron changes its momentum. In that sense, also the phonon interaction 
leads to an infinite hierarchy formally equivalent to the many-body hierarchy problem discussed above. For example, 

^{Dp,p^a\ fe'^A,/c-p) = {Dp.p^a\ kO-x.k-p) — {Dp,pj}{a\ k^>^-k-p) describes correlations between carriers and phonons. 
A similar treatment of the cuiantized light-matter interaction has led to quantum-optical effects resulting, e.g., in 
squeezing and entanglement .Ejo'lj The singlet-doublet truncation with a controlled scattering treatment of the 
triplets leads to a closed set of equations which includes the dominant correlations between lattice vibrations and 
carriers. 



Coulomb Interaction 



Since the Heisenberg equations of motion are linear with respect to the different parts of the total Hamiltonian (||), 
we examine these different parts separately. For the exciton correlations under the influence of Coulomb interaction, 
we obtain 



d 



f- ^ q,k' ,k 

in—c 



~h 



~h 

^k-q 



q,k ,k 



+ Vk^k'-q [(1 - fk) (1 - fk-q) fk'+qfk' fkfk-q (l ^ fk'+q) (l fk')] 
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k' —q 



{fk-q Ik') 



^q+k'-l,Lk 
"e+X 



{Ik fk'+q) 



^q—k+l,k',l 



q,k' .1 



q.l,k 

'x 



q+fc —l,k — q-\-l 



[fk'+q~ft\Y.'^l-kd 



q-~k-\-Lk' A 



(16) 



which is coupled both to carrier densities and to electron and hole correlations c 



q,k ,k 



via the terms ■^'^'^ 



-e(h) "-x The dynamics of these correlations is described by similar equations given in App. @. 
The interpretation of Eq. ( [l6| ) is straightforward; the first hne gives the kinetic evolution of the four-point operator 



Q,fc ,k , q.k ,k 



with the rcnormalized energies et''*'* = e 



e(h) 



J2k' Vk^k'fl, 



The second line contains the factorized source term 



which initiates the creation of excitonic correlations as soon as electrons and holes are present. This singlet source 
is altered by the direct influence of the correlations in the third line. The remaining six Coulomb sums describe the 
different possibilities how two out of the four fermions in A(aJa|^Oca„) can interact via the Coulomb interaction. These 
sums lead to the possibility to form bound excitons. The major contributions originate from the first two of the six 
sums which are multiplied by a phase space filling factor instead of a density difference and are thus appreciable even 
for low densities. All other sums vanish for low density but become important when the density is increased. They 
are a consequence of the indistinguishability of electrons and holes and correspond to Coulomb interaction betweea, 
carriers formally "belonging" to two different excitons. Equations investigating exciton formation in the x^^^-liniitlii'ta 
do neither contain these additional /-dependent contributions nor the singlet source. 

In order to obtain a closed set of equations for the incoherent dynamics of the pure carrier system, we also have to 
derive the equation of motion for electron and hole densities. They are given by 



Im 



E^. 



q^e 



k' ,q 



k' ,q 



-Im 



EVf +Y.^k-q^^'C^ 

k' .q k' ,q 



q^k.k' 



(17) 



(18) 



and describe the influence of the correlations on the carrier distributions. Since all incoherent four-point correla- 
tions are treated exactly, this approach includes the microscopic carrier-carrier scattering beyond the second Born 
approach. El 



D. Phonon Interaction 



As a consequence of the Coulomb interaction, excitonic correlations can build up since the carriers of opposite 
charge are attracted towards each other. However, two particles cannot be bound without the assistance of a third 
object due to energy and momentum conservation. In a many-body system, any additional carrier can be this third 
object and compensate the energy gained by the bound-state formation. This process inevitably leads to heating of 
the remaining carrier system such that further formation gets more and more improbable. By a coupling to phonons 
the excess energy can be directed out of the many-body carrier system into the phonon reservoir. For simplicity, we 
assume that the phonon system is basically unperturbed by the carrier dynamics and solve the phonon dynamics using 
the steady state Markov approximation with the assumption of thermal occupation of the different phonon states. 
Working out the phonon contributions to the dynamic equations, we obtain in the singlet-doublet approximation 

= lY.^m[UlJ , (19) 

Hp p 

= -|Ei^^^[ny . (20) 

Hp p 

= ^k,k-q-k' {fk' ^ fk-q) ^ ^k',k'+q-k {fk ^ fk'+q) : 

Hp, SB 




aj-J k 



■t ^ q,k',k 
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' dt ^ 



■ * d q k' ,k 

in—cV 



n 



dt 



ffp,SD 



Hp.SD 



k,k—q—k 



• {fk-q ^ fk') + ^k',k'+q-k {fk'+q ~ fk) 



+ ^k',-q {fk fk-q) + ^k,q {fk' ^ fk'+q) i 
^k,k-q-k' {fk' ~ fk-q) + nfe/^fc'+q-fe {fk ^ fk'+q) 
+ ^k',-q {fk-q fk) + {fk'+q ~ fk') 



(21) 



(22) 



(23) 



where we have defined the term 



k—p 



nn 



p± 



p± 



■ y ^ \^p,p± I 



{5(ef+p - ef + 7jf^p,pi ) - gief+p - ef - ?if^p,pi ) } , 



|2 k—p~l,l,k 



'i+P 



- nn 



P:Pi)}i 



(24) 



and a similar expression for with (7(e) = 'KS{e) + iP{-^)- The derivation of Eqs. (|l9|)-(|24|) is presented in App. 
After being inserted into Eqs. ([1^) and (^), the first four hnes of Eq. ( p^ ) have a straight forward interpretation as 
scattering rates. They are the typical scattering equations with a balance between scattering into (lines one and two) 
and scattering out of (lines three and four) electron state k. In both cases, phonon absorption and emission processes 
are possible. Since emission can happen even without any phonons present, the respective terms are proportional to 
Np^^ + 1. The last two lines are due to phonon coupling to higher order correlations. Since both the carrier and the 
correlation dynamics in Eqs. ([l9|)-(p3|) are determined by similar H-terms, the corresponding phonon effects result 
from the cooling of the carrier plasma. 

The actual exciton formation dynamics is described in the six-point (triplet) level. The major contribution of 
these triplet correlations can be evaluated at the scattering levelEj in analogy to the second-Born approximation 
four-point correlations which has successfully been used to describe microscopic Coulomb and phonon scattering 
Physically, this approach provides scattering of excitons with a third object excluding the formation of bound three- 
particle clusters. As discussed in App. the resulting equations can be written as 



•fc ^ q,k',k 

dt ^ 



Hp,T 



p 

-Efe-p-(^fe%,fe-p)*)4'''" 



E(7^,- 



iik' 



k' -\-q,k' ■ 



-pT 



<l:P-k 



-E((^t..p-«)* + (^^^'+.,.-p)*)4''' 
p 

+ EW',p-,+7.Vp) (4'-^-'='+')* 
p 

^ E {^'k',p-q ^ ilk-q,p-q)*) '^X 
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j).k' ,k — q-\-p 

'x 



(25) 



where 7 and 7 are defined according to 



7L 



p± 

= (~z) ^ |G,,,, p [(7v™^ + /^')5(etp - £^ + nn^^.j 
+ (ivP«^ + 1 - f^)g{stp -el- nn,^,j] , 

= (-Z) 5] |Gp,p, p [(7V™^ + - + Tif^p.pJ 

Pi 

+ + 1 - /fe-p)5(4 - 4-p - 7ir!p,pj] , 

= (-z) ^ |Gp,p, p [(TVPH^ + 1 - ft-M^k^P ' + ?if^p.p J 
+ + /tp)5(etp - 4 - Tt^^p^pJ] ■ 



(26) 



(27) 



(28) 



(29) 



Similar equations for electron and hole correlations are also given in App. 

These six-point phonon scattering contributions lead to microscopic dephasing of correlations and we find the 
conservation law 



k,k',q 



q,k\k 
A 



(30) 



Hp 



for the complex valued correlations c'^^ such that the phonon scattering has a diffusive character. This property 
allows formation as well as equilibra tion of correlations. 

The set of equations (|l6|)-(|l|), ( |Ai1 ), and (|a|), extended by the phonon contributions, Eqs. (|9l)-(|2|), (^ ) 
and (B8), provides a closed system of coupled differential equations which can be solved under different initial condi- 
tions in order to study exciton formation for different carrier densities and lattice temperatures. 



III. EXCITON FORMATION 



Because of the high dimensionality of the summations in our hierarchy of coupled equations, the currently available 
computer resources limit us to evaluate the equations for a one-dimensional model system. However, we choose 
parameters to be close to the standard GaAs parameters used for quantum wells. The effective width of our quantum 
wire is determined such that the exciton binding energy lies roughly llmeV below the unrenormalized bandgap 
which is also the case in typical 8nm quantum wells. The corresponding 3D-exciton Bohr radius is ao = 12.4 nm. 
To incorporate Coulomb screening effects to the four-point terms, we use a statically screened Coulomb potential 
obtained from the Lindhard formulae Microscopically, the justification follows from the coupling to the six-point 
terms. For the carrier densities, the coupling to the doublet level provides a microscopic description of scattering and 
screening such that the unscreened matrix element is used in Eq. ( p^ and (|l^). In all calculations, we start from an 
incoherent electron-hole plasma with vanishing correlated doublets. 



A. Direct Evidence of Exciton Formation 



The formation of-aiw significant amount of bound excitons can be directly observed via the electron-hole pair- 
correlation functioncj'Eil 

g^^r) EE ^{^l{r)¥,imhme{r)}, (31) 

which is normalized with respect to the single-particle carrier density n and measures the conditional probability to 
find an electron at position r while the hole is located at r = 0. The pair-correlation function can be separated into 
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a singlet and a correlated doublet part 



^ ' k,k',p,p' 



nC 



^ ' k,k',q 



-i(k~k'-q)r q,k',k 



(32) 

(33) 
(34) 



where Pk = (aj, f.ac,k) is the coherent microscopic interband transition amplitude which provides the dominant contri- 
bution to the pair-correlation function after a coherent excitation with a classical field.E3 As soon as this polarization 
dephases, however, the factorized contribution becomes constant and the only r-dependence is provided by the corre- 
lated part Ag'^'*(r). This quantity provides an intuitive measure of when and how fast exciton formation and ionization 
takes place in incoherent configurations. 




FIG. 1: Pair-correlation function Ag''^{r)/{rfn^) normalized to the constant Hartree-Fock value for a lattice temperature of 
T = 10 K and a carrier density of n = 2 x 10* cm~^ at different times. For comparison, the wave function of the lowest bound 
exciton is given as a shaded area. The dashed line is the corresponding result of a second-Born computation. The inset shows 
a magnification of the tails of the same curves. 



Figure^ shows computed examples of Ag'^^{r) as a function of electron- hole distance r for a low carrier density of 
10 cm . At early times after the completely uncorrelated initialization of the system, the carriers react to 
the Coulomb attraction and the probability of finding electrons and holes close to each other increases. However, the 
correlated Ag'^^ has clearly negative parts. At certain distances, the total pair correlation g'''^(r) is therefore below its 
factorized value. We interpret this behavior as a fast rearrangement of the electron-hole plasma where the probability 
to find electrons in the close vicinity of holes is increased at the expense of |tiie probability at around 2 to 4 ao. Since 
the area under Ag'^'* provides a measure of how many excitons have formedcj, the early time dynamics does not show 
a real formation of excitons out of the plasma. Nevertheless, at 480 ps, the pair correlation has assumed the shape of 
an exciton relative-motion wave function (given as a shaded area) and grows almost linearly in time. 

For comparison, we also show the second Born result in Fig. [I] as a dashed line. In this approximation, the Coulomb 
sums of Eqs. (|l^) and the triplet phonon scattering, Eq. (|2|), are not included such that the doublets can describe 
only carrier-carrier scattering but not bound excitons. We see that the resulting pair-correlation function has similar 
characteristics as the early time Ag'^'^ in the full computation. In particular, Ag^^ drops to negative values before it 
approaches zero at large distances. Since no exciton populations are included in the second-Born computation, the 
resulting shape of the pair correlation has to be interpreted as a rearrangement within the plasma. 
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B. Phonon Induced Energy Transfer in Exciton Formation 

Since a bound exciton at rest has a lower energy than a free electron- hole pair, we study the energy transfer during 
the exciton formation in the many-body system. For carriers alone, the system energy is 

{Hkin) + {Hc) = J2^k{a[^k^\.k) + ^ ^ii^l,k<^\',k'^>^'^k'+qax,k-q) 

X,k A, A' 

fc,fc' ,q^0 

+A{H^c) + + mc) (35) 

with the kinetic energies (H^^^) — J^k^kfk' ^^'^ incoherent singlet (factorized) parts to the Coulomb energy 
{H^)s = -h^Ek^'^kVk-k'f^J^, and the correlated contributions A{H^) = i Efe,fc',g^o ^cf and A(iJg''} = 
-EkM.,^oVA''''' Using Eqs. we find 

^ [(ifkin) + {Hc)]h,,„+Hc = (36) 

which shows that the energy is conserved within the isolated carrier system. As a result, exciton formation is not 
efficient for a plain carrier system. When also phonons are included, we find 

QliHtot) = ^ [(i?ki„) + (i^c) + (^fphon) + {Hp)] = , (37) 

i.e., the truncation via the cluster expansion scheme fully conserves the total energy of the system. When the phonons 
are included, the carrier system can cool down because part of the energy gained via exciton formation is transferred 
to lattice vibrations. When the phonons are treated as a reservoir, the phonon bath acts as a sink and the energy 
flux directed to the phonon system is absorbed by the reservoir. However, when we replace the microscopic phonon 
scattering, Eq. (p5[), by a simplified constant dephasing rate 7, we find 

^(i?)tot = -7A(i?c) >0 (38) 

which implies an unphysical heating of the system. Therefore, systematic exciton formation studies require a micro- 
scopic description of phonon scattering. 

The dynamic evolution of the carrier energy is shown in Fig. ^ by using either constant dephasing or microscopic 
phonon scattering. With constant 7, the attractive Coulomb energy very quickly reaches its steady state value such 
that significant pair formation is not observed. Even a small constant 7 of 50 /leV leads to a linear heating of the 
carrier distributions, i.e., to an increase of the total energy. With the microscopic phonon scattering and a lattice 
temperature of 10 K, however, true formation of excitons is possible as indicated by the continuous decrease of the 
Coulomb energy. Simultaneously, the total energy of the carrier system decreases since part of the energy is lost to 
the phonon bath. The slope of decrease of the attractive Coulomb energy can in principle be used as a measure of 
how fast the exciton formation takes place. 

To speed up the numerics, we have computed Fig. |^ by following the dynamics upto 40 ps with a phonon matrix- 
element enhanced by a factor of 30; this corresponds to 30 x 40ps = 1.2 ns of formation dynamics. To verify the 
validity of this analysis we show in Fig. ^ a comparison for three enhancement factors of 10, 15, and 30. After rescaling 
the time axis with the respective factor, the behavior of the correlated Coulomb energy is indeed independent of the 
enhancement factor. In particular, the rate of change is very similar for all three cases which justifies the use of the 
artificial enhancement factors in the numerical investigations. 



C. Formation of Specific Excitons 

In the equation for excitonic correlations Eq. (|I^), two of the six Coulomb sums are usually dominant; the terms 
multiplied by the phase space filling factor 1 — /"^ — provide a large contribution for all densities. Therefore, the 
restriction to these two dominant sums is often a very good approximation to the full result. The corresponding 
"main sum approximation" of Eq. (|l6|) is 

r) 

Afj_„<]:k ^ (?e ~h _?e,~h\q.k,k 

^'^ Q^'^X —\^k'+q ^k-q ^k + ^k')^X 
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FIG. 2: Attractive Coulomb energy (a) and total energy (b) per particle comparing a computation with constant dephasing 
approximation (dashed line) with the full result including microscopic phonon scattering (solid line) for the same parameters 
as in Fig. ]l. Energies are given in multiples of the 3D exciton binding energy Eb = 4.2 mcV. 
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FIG. 3: Correlated Coulomb energy A{Hq ) per particle for three different computations with different enhancement factors. 
For each curve, the time axis has been rescaled with the respective enhancement factor. The carrier density is n = 1 x 10* cm^^ 
and the lattice temperature is T = 10 K. 



+Vk-,,^, [(1 - /|)(1 - ft,)fk'+Jk' - (1 - - fk')fkf'k-,] 

+ [1 - fk - ft,] E ^-^4'' ' " [1 - fk'+, - fk'] E Vi-k'ct' ■ (39) 
I I 

With this restriction, the excitonic correlations are not coupled to carrier-carrier correlations anymore and we obtain 
a closed subsystem of equations for each center-of-mass momentum q. 

Since the carrier distributions typically vary slowly, we try to find an adiabatic solution for Eq. (|39|). To do so, we 
use a generalized exciton basis, which is defined by the equations 

ek,,^:,,ik) - (1 - - /t,.) E Vk'-k^:,,{k') = E,,,^:jk) , (4o) 

k' 

(<,(fc))*eM-E(<<?(^'))*^^-'-fe(l-/fc'+r -/fe'-,'^) = (41) 

k' 
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where we have defined ek,q = ^t+q^ +^k-q'^ l'^^'^^ ~ '?"^e(h)/(™e + ™/i)- The main difference to the usual Wannier 
basis is that for non-zero densities the phase-space fiUing factor and a screened Coulomb potential enter the effective 
Hamiltonian. Therefore, the eigenvalue problem Eq. ( |40| ) is not Hermitian and one obtains left and right handed 

eigenfunctions (^y^q . This exciton basis fulfills a generalized orthogonality and completeness relation 



k 



k,k' 



(42) 
(43) 



The center-of-mass momentum q enters Eq. (40) via the phase-space filling factor such that the relative motion of an 
exciton depends also on its center-of-mass momentum. 

By means of the generalized Wannier functions, one can introduce exciton annihilation operators 



k 

vl_qHCk + q^ = ^(K^qik) Xy^q 



(44) 

(45) 



and the Hermitian conjugate creation operators. Equations ( ^ ) and (45) can directly be applied to the correlations 
in order to transform them according to 



q,k'-q .k+q" _ 



J2{^:jk)yq^:,jk')A{xi^Xy,^q), 



A{Xl^Xy,,q) = ^<,(fc)(<^L,,(fc')) 



'W* ^q.k'-q'\k+q'' 



(46) 
(47) 



k,k' 



This expansion cannot only be used to solve Eq. ( |39| ) but also to compute expectation values in the exciton basis from 
our numerical computations performed in the fc-basis. 
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FIG. 4: Formation dynamics of q- integrated population correlations, Eq. (^sj), out of an electron-hole plasma with n = 
2 X lO^cm"^ at T = 60 K. The lattice temperature was 10 K (solid line), 20 K (dotted line), 30 K (dashed line), and 40 K 
(dash-dotted line), respectively. Thick lines indicate A{A'^is) whereas thin lines refer to A(A'2p). 

Using Eq. (|4^), we can study the formation of specific excitons as shown in Fig. |^. There, we have computed the 
dynamics for different lattice temperatures, always starting with initial electron and hole distributions at a temperature 
of T = 60 K. With the help of Eq. (E^), we compute the g- integrated exciton population correlations 



A(iV.) = ^ A(iV.,,) = J2 MXlqX^^q) 



(48) 
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normalized to the total number of electrons in the system. ICJiese populations as a function of time are shown for the 
two lowest bound excitons denoted Is and 2p, respectivelyEJ. Figure ^ clearly demonstrates how the formation rate 
of Is excitons quickly drops for elevated lattice temperatures. For temperatures above 40 K, we do not expect any 
significant formation. 

In order to continue our analytical derivation, we note that in the incoherent regime and typically change 
slowly such that the exciton wave functions can be assumed to be quasi stationary. With the help of this adiabatic 
approximation, Eq. (p59) can be written as 



(49) 



k,k' 



with the factorized (singlet) source term 



(-'- fk+q" fk-q'^)fk'+q''fk'-q'^ 

~ (1 ^ fk'+q" ^ fk'-q'^)fk+q<^fk-ql^ 



(50) 



The properties of the left-handed exciton basis, Eq. (^), can be used to simplify the last term in Eq. ( |4S| ) according 



to 



E<.W(<^^',.(^'))*^''''''' = {E.--E,){Xl^X,.^,), 



(51) 



k,k' 



where the factorized plasma part of the two-particle correlation is given by 

{Xl,Xu',q)s = ('^^'..(^))* fk+q'fk-q- 



(52) 



In the exciton basis, the full equation is therefore 

d 



ih-A{Xl^^X,,,g) = (E,, - E,)A{Xl^X,y^,) + {E,, - E,){Xl^X,, 



-A{Xl^X,,.,) 



(53) 



T-scatt 



where we have included the scattering contributions from triplet correlations. 

The simple form of Eq. (^ ) reveals how the exciton formation proceeds. If we start the calculations assuming 
initially vanishing correlations, they start to build up because the source term {Xl ^X,^/ ^q)^ is non- vanishing as 
soon as carriers are present in the system. However, that source term does not create diagonal correlations such that 
excitonic populations A{N^,q) stay zero without the six-point scattering contribution. Hence, the formation of exciton 
populations in the incoherent regime has only one possible Coulombic channel where off-diagonal transition correlations 
A{Xl qXt,>^^,q) are created first. Only after off-diagonal transitions have been created, populations can be formed 
via the triplet scattering. Therefore, it is crucial to include both excitonic transition and population correlations 
in the formation analysis, as is done automatically in the full fc-basis treatment used in all our computations. If 
one uses excitonic models which are restricted to the diagonal populations, the formation channel via the transition 
correlations is completely omitted. 

The efficiency and the character of the indirect exciton formation channel are analyzed in Fig. ^ which investigates 
the build-up of the integrated exciton population correlation A(A^is). Here, we compare a run where all off-diagonal 
transition correlations A{X'lX^i^i^) are set to zero during the computation with a reference run without reset. For 
the reset, all correlations are transformed from the fc-basis into the exciton basis via Eq. ( ^7|) and back to the fc-basis 
with the help of Eq. (^). In the back transformation, we omit all but the diagonal correlations. No phonon- 
enhancement factor was used in this computation because we want to investigate Coulomb and phonon dynamics 
for equally short time scales. Exactly at the moment where the off-diagonal correlations are reset, the formation 
dynamics stops and no more exciton correlations can build up. On the contrary, the amount of Is populations is 
even slighly decreased since the phonon scattering tends to transfer some of the diagonal populations back to off- 
diagonal transition correlations. Already around 3ps after the reset, however, the off-diagonal transition correlations 
are fully recovered. After some transient dynamics, the formation continues with the previous formation rate. Also 
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FIG. 5: Total exciton population correlation A(Afis,5) for two computations with (dotted) and without (dashed) reset of 
off-diagonal transition correlations for n = 2 x 10'' cm~^ and T — 10 K. The thin solid line is given as guide to the eye. 



in other quantities, as for example in the correlated Coulomb energy, this fast recovery of the transition correlations 
is confirmed. This behavior indicates that off-diagonal transition correlations build up extremely quickly compared 
to the formation rate of diagonal populations. It is these off-diagonal transition correlations which cause the plasma 
like shape of the early-time pair-correlation function. 

The importance of the off-diagonal exciton transitions can also be seen in the equation of motion of the electron 
density. By expressing the right hand side of Eq. (^^ with the help of Eq. (^6|), we obtain 



d_ 

dt 



X— corrs 



(54) 



Thus, diagonal exciton populations do not change the carrier densities. Excitonic transitions are required in order to 
describe the full dynamics and the correct heating of the carrier plasma due to exciton formation. 

D. Transition Correlations 

In order to understand the general nature of the off-diagonal transition correlations, it is sufficient to investigate 
Eq. ( p3| ) where the six-point scattering contributions are approximated by a constant dephasing rate, i.e.. 



J T-scatt 

Since this term leads to a simple decay, Eq. (|5^) has a steady-state solution 



For sufficiently small 7, we find 



Eu' — Ed 



steady —state 



Eu - 

A{Xl^X,,^,) = {6,y - l){Xl^X,,^,)s . 



q/S ■ 



(55) 



(56) 



(57) 



For v^v' , the off-diagonal correlations A(X^ ^X^i ^q) — —{X^ ^X^i ^q)^ fully cancel with the singlet (two-point) contri- 
bution when the total {X^ ^X^i ^q) is evaluated. The diagonal part {Xl ^X,y^q) , however, is determined by the singlet 
contribution alone. In other words, we find 



{^l,q^y'.q) ~ 5y,y{Xl^^Xi,^q)s = 5y,y ^ Wu,q{^)? fk+q" fk-qf- ' 



(58) 
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FIG. 6: Exciton correlations (JfJX^) as function of energy Ei, after 720 ps of computation for the same parameters as in Fig. |^. 
The inset magnifies the region for energies higher than Eg — I.SEb- While off-diagonal correlations are largely reduced in the 
full computation, the diagonal exciton population increases in time due to genuine exciton formation. 



where the solution is expressed with the help of Eq. ( p^ ) . 

Clearly, the ofF-diagonal transition correlations build up to compensate the off-diagonal two-point contribu- 
tions whereas diagonal populations {Xl gX^^q) have a non-vanishing and observable fermionic plasma contribution 
{Xl qX^i ^q)s. These observations are valid also for the more general case where the six-point phonon scattering is 

included microscopically. Figure ^ shows the absolute value of the total exciton correlations \{xl^ q=o^iy.q=o)\ a 
function of the exciton energy E^, for such a general calculation with microscopic phonon scattering. We compare the 
results of the full calculation with those obtained if we truncate the equations at the singlet level where no excitonic 
correlations are included. 

The marker at the lowest energy gives the diagonal Is exciton population whereas the other markers indicate the 
magnitude of the respective off-diagonal transitions. We see, that a calculation at the singlet level predicts that the 
contributions of the off-diagonal transitions are quite appreciable (open squares) whereas the full calculation yields 
results that are more than an order of magnitude smaller (full circles). Hence, we find a high degree of cancellation 
between the singlet and the higher-order results as predicted by Eq. (^8|). Due to the very fast build-up of off- 
diagonal correlations compared to the slow phonon scattering, the steady-state result Eq. ( |57| ) for off-diagonal exciton 
transitions is still a good approximation for the full calculation. The main change to the simplified analysis leading to 
Eq. (|5^) is that under true formation conditions the diagonal exciton populations can grow and exceed their singlet 
contribution. Due to the formation dynamics, this diagonal population may increase in time and the difference to the 
plasma level determines the amount of true populations. This difference can also be seen in Fig. ^. 

Even though the final correlations {X^ qX,^i ^g) are to a very good approximation diagonal in the exciton basis, the 
distinction between factorized and correlated part is very useful because the consistent theory shows that: i) only 
the truly correlated part of the off-diagonal exciton correlations leads to a significant contribution to the formation of 
excitons out of an incoherent plasma, ii) only the off-diagonal correlations influence the time evolution of the carrier 
distributions, and iii) even under good formation conditions, the correlated contribution to the exciton population is 
still of the same order of magnitude as its factorized counterpart. Both contributions certainly influence experiments 
and must be carefully distinguished. FctLjexample, the excitonic photoluminescence does not require the presence of 
true excitonic populations A{Xl gXiy_q)S3 



IV. SUMMARY 



We have presented a general many-body formalism to describe the dynamics of charge carriers and two-particle 
correlations within semiconductor heterostructures. The Coulomb interaction between electrons and holes and the 
coupling of carriers and correlations to a phonon reservoir has been described microscopically. We thus are able 
to compute exciton formation times and exciton-exciton correlations without neglecting the underlying fermionic 
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properties. We have shown that the formation is an intricate process where Coulomb correlations rapidly build up 
on a picosecond time scale while phonon dynamics leads to true exciton formation on a slow nanosecond time scale. 
An adiabatic approximative solution has been obtained which provides an intuitive interpretation of the microscopic 
results. 
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FIG. 7; Contour plot of exciton formation rates for an initial carrier temperature of 60 K and varying lattice temperatures 
and carrier densities. The exciton formation rates are taken at 1.2 ns of evolution and normalized by the carrier density. The 
contour lines thus determine which percentage of charge carriers are bound to Is-excitons per nanosecond. 

The basic results are summarized in Fig. ^ There, we show the exciton formation rate 1.2 ns after initialization of 
the computations. At this time, the formation proceeds almost linearly in time, c.f. Fig. ^, and the formation rate gives 
a good measure of the efficiency of the exciton formation. The contour lines show how many percent of electron- hole 
pairs are bound to Is-excitons per nanosecond. Exciton formation is only efficient for cold lattice temperatures below 
30 K and rather low densities. For too low densities, however, the rate of formation drops simply because it becomes 
increasingly improbable that electrons and holes find scattering partners necessary for the formation process. 

Generally, we can interpret the formation of excitons as a build-up of correlations which does not change the number 
of electrons and holes in the system. An interpretation of excitons as atom-like entities only makes sense under very 
restricted conditions of a dilute gas where a small density of excitons interacts only very weakly with the remaining 
carriers. Under those conditions, our theory gives reasonable exciton numbers. Basically, the total number of excitons 
formed is controlled by the lifetime of the electron-hole excitations. In many realistic systems, this is in the order of 
ns because of non-radiative recombination at imperfections and/or impurities. 
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APPENDIX A: EQUATIONS FOR ELECTRON AND HOLE CORRELATIONS 



In Section II C we have only presented the equation for the excitonic correlations 
for electron-electron and hole-hole correlations are 



q,k' ,k 



^k-q + ^k'+q 



^q.k' ,k 



q [fk-qfk'+q (1 ^ fk) (1 ^ fk') 



-Vk- 

-Vq [fLqfk'+q (1 fk) (1 fk') fkfk' (l " fk-q) (l fk'+q)] 



The remaining equations 



fkfk' (1 - fk-q) (1 - fk'+q)] 
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-V„ 



{fk'+q fk') C^eA 



q-ld,k 
X 



{fk fk-q) 



^q+k' -l,l,k' 



-Vk-k' 



{fk'+q fk) 



k — q—l,l.k 



/-/■e re g+fc'-/,i,/c 



[ft - /L,] E ^-^-c?''' ' + [AVg - fk] E ^/-^^cr'+'''''' 
[ft'+q - /fc'] E ^-fe'c^''' - [fi - fk-q] E yi-k'cr"'- 



-l.Lk 



(Al) 



2h I q,k' ,k 



+ l^.-.'-g [(1 - ftq) (1 - fk' + q) fkfk' - (1 - (1 - fk') fk-qfk'+q] 

-Vq [(1 - ftq) (1 - Z^'+g) fkf'k' - (1 - (1 - fk') fk-qfk' + q] 



+ V„ 



{fk' fk'+q) E 



q—k+l.k.l 
'h+X 



{fk-q ^ fk) E 



J — q — k',k'.l 
"^h+X 



+ Vk-k'-q 



{fk'+q fk) E 



q—k+l.k'^l 
"^h+X 



{fk' fk-q) E "^e+Jf 



l — q — k',k,l 



[1 - /tg - /^Vg] E ^-.^5;'' ' + [1 - - fH'] E ^-.^'^ 



[/' - fk-q] E ^-^4''' ' - [/^'+g - A'] E ^'-'^^^ 



q — k-\-l,k' ,1 

[fli'+q - fl^ E "^^-^'^t" + [ft' ~ ftq] E Vl-k'Cr'-'''t 



(A2) 



APPENDIX B: PHONON INTERACTION 

The general phonon interaction Hamiltonian Eq. (Q) is the starting point to compute the operator equation of 
motion for a two-point operator 



ih 



d 



Hp 



' Q^-^\,k^^,k+q 

where the collective phonon operator 

Gp — E] ^p.Pi y^v, 



= E^P («A,A:«A,/c+g-p - al,fc+p"A,/c+g ) , 



TV 

P± ^ ^-'P,pi_ 



(Bl) 



(B2) 



pi- 



has been defined. With the help of this operator equation, we can now generate all other equations of motion. The 
simplest case is the equation for carrier densities. For electrons, we obtain 



Hp 



E ^{Gp"-l,k"'c,k-p) - h-c. 



(B3) 



In order to compute the right-hand side of Eq. (B3), we have to establish the equations of motion for the phonon 
assisted terms. In the case of the phonon assisted electron density, this gives 



d 



i^-Q^MDp,p^al,^ac,k-p) 



— Gp^p^fk-p{l — fk) 
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-APp,p,(7_p)(/,^-/,%) 



k—p—l.l,k 



(B4) 



where we have aheady left out all terms proportional to the coherent transition amplitudes Pk = {v],Ck)- In cases 
with coherent excitations, those terms can easily be included. But since in this case the Coulomb scattering is usually 
dominant, th e m ain application of the microscopic phonon scattering lies in the incoherent regime. In order to 
compute Eq. (|B3|) , we solve Eq. (B4) in Markov approximation. In principle, the result of a Markov approximation 
is dependent on the basis which we use . In general, however, the difference is not crucial if the resulting terms are 
summed over as is the case in Eq. (B3). For the computations, we have therefore included only terms proportional 
to the square of the phonon matrix elements. No Coulomb or light matter interaction was included in the equation 
of phonon assisted densities. The Markov approximation was done in the one-particle basis of Bloch electrons and 
holes. We numerically confirmed that our result is independent of the precise choice of the basis. Additionally, we 
assume a reservoir of phonons, i.e., we neglect coherent phonons and set {Dp^p^Dp,pj) = (exp(-Ej,,p^ /(fcT)) — 1) ^ 
to the Bose-Einstein distribution at the corresponding energy Ep p_^ given by the phonon dispersion. After those 
approximations, we obtain Eqs. (|l9|) and (|20| ) for the carrier densities. 

We have seen in Sec. II D that Eq. ( |B4D and its counterpart for the phonon assisted hole densities suffice to 
compute the doublet contributions to the correlation equations. In addition, phonon assisted correlations of the form 
were needed in order to provide a true dephasing mechanism for excitons and carrier-carrier correlations 
on the triplet level. More specifically, the phonon contribution to the equation of motion for excitonic correlations is 



d 

iT^QlHalkalf^,ac,k'+qav,k-q) 



-^{Gl'+q-k0^ik^c,k'+q) {ftq fk') " ^{Gk' +q-kalk' ,k-q) [ft " fk'+q) 

- E A{glal^,^al yac,k'+qav,k^q+p) ~ ^ ^{Qlal^k-pO'l,k''^c,k'+qav,k-q) 



H ^^'3lal^kal,k'0'c,k'+q+pav,k-q) - ^{Gl 



,k' -p^':,^' +qO'V ,k-q) 



(B5) 



and similar equations exis t also for the carrier-carrier correlations. The first row involves the phonon assisted densities 
solved already in Eq. (B4). The remaining phonon assisted correlations are solved via 



d 



Gp,p^ (^p™^ + fk-q-p) /^{al kalj^.ac.k'+qavM-q) 

+ Gp,p± {Np^^_ + 1 ^ fk'+q) ^{O'c.k'^vM'^'C,^' +q+pav,k-q-p) 

-Gp.P^ (A^p™ + 1 - fb) A(a^.al;,fc,_pac,fc'+,a.,fe-,-p) 

-Gp,p^ (^™^ + fk) M(i,k-p(^Vk'"-c,k'+q<^vM-q-p) 



(B6) 



which together with similar terms for phonon-assisted carrier-carrier correlations can be used to build up the complete 
phonon interaction contributions to the correlation equations. Again, we have kept only incoherent terms which enter 
the correlation equations proportional to the square of the phonon matrix element. Equation (B6) is solved in Markov 
approximation assuming a phonon bath characterized by the lattice temperature. This way, we obtain Eq. ( p5| ) and 
the corresponding equations 



Hp,T 



p 



E (7fc',fc'-p 



\ ' k +q.k —p) 



rl'P'k 



{{ll^q^p^qY + ill'+q^q^pY) ^ + E ^p-q + 7.%-p) (cf +^ ) * 



P 

E 

p 

ipik' ,p-q 

P 



\ lk — q,p—q) ) / , \ lk,q—p \ lk'+q,q~pj ) 

P 



p,k' ,k-q+p 



(B7) 
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and 



■t ^ q,k',k 



<x,h \*\ ^q^k'.k 



p 

~ E {'^k,k-p ~ (7fc-<j,A;-p)*) (^h ~ E (7fc',fc'-p ~ (7fe'+q,fe'-p)*) C^'''' 

P P 

- E H-,,,-qr + (7^+,,,-p)*) cf'^' + y: {i'k',,-q + 7t-p) {c'-'-'"''^r 
p p 

- E (7^,P-. - (7t,.p-,)*) c^/-'+^-^-^- - E (7L-P - (7.V,,,-p)*) c^/-' 
p p 

where the terms 7 and 7 are identical to those defined in Eqs. (]26|)-(p9|). 



(B8) 
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